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ABSTRACT 

It  is  shown  that  a  bounded  biinfinite  banded  totally  positive  matrix 
A  is  boundedly  invertible  iff  there  is  one  and  only  one  bounded  sequence 
mapped  by  A  to  the  sequence  ((-)1).  The  argument  shows  that  such  a 
matrix  has  a  main  diagonal,  i.e.,  the  inverse  of  A  is  the  bounded  point- 


wise  limit  of  inverses  of  finite  sections  of  A  principal  with  respect  to 
a  particular  diagonal,  hence  {(-)1+^A  ^(i,j))  or  its  negative  is  again 
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SIGNIFICANCE  AND  EXPLANATION 
Spline  approximation  is  often  most  effective  when  the  breakpoint  (knot) 
sequence  can  be  chosen  suitably  non-uniform*  At  the  same  time,  the  standard 
approximation  schemes  (such  as  least-squares  approximation,  or  interpolation 
at  suitable  interpolation  points  by  splines)  are  so  far  only  known  to  be 
usable  and  bounded  as  long  as  the  breakpoint  sequence  is  almost  uniform. 

The  problem  of  showing  existence  and  uniqueness  of  bounded  spline  approxi- 
mants  to  bounded  data  boils  down  to  showing  invertibility  of  a  certain 
infinite  matrix  A.  The  distinguished  features  of  this  matrix  are  its 
bandedness  and  its  total  positivity,  i.e.,  all  minors  of  A  are  nonnegative. 
In  this  paper  we  show  that  if  there  is  exactly  one  bounded  sequence  mapped 
by  a  bi infinite  totally  positive  banded  matrix  A  to  the  particular 
sequence  (• • • ,-l, 1,-1, 1,-1,* ••) ,  then  every  bounded  sequence  is  contained 
in  the  range  of  A.  In  spline  terms,  this  result  says,  for  example,  that 
any  bounded  data  sequence  can  be  interpolated,  and  in  exactly  one  way,  with 
a  bounded  spline  (with  a  given  knot  sequence,  at  a  given  interpolation  point 
sequence)  provided  that  the  periodic  data  ^+1,-1^  can  be  interpolated,  and 
in  exactly  one  way,  by  a  bounded  spline  from  that  class.  Further,  our  argu¬ 
ments  show  that  such  an  interpolating  spline  can  be  constructed  as  the  limit 


of  splines  which  satisfy  finitely  many  of  the  given  interpolation  conditions. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


THE  INVERSE  OF  A  TOTALLY  POSITIVE  BIINFINITE  BAND  MATRIX 
Carl  de  Boor 

0 • Introduction •  This  paper  is  a  further  step  in  a  continuing  effort  to  understand  cer¬ 
tain  linear  spline  approximation  schemes*  Convergence  of  such  processes  is  intimately  tied 
to  their  stability,  i*e*,  to  their  boundedness,  as  maps  on  C  ,  say*  Use  of  the  B-spline 
basis  shows  this  question  to  be  equivalent  to  bounding  the  Inverse  of  certain  totally  posi¬ 
tive  band  matrices*  The  calculation  of  bounds  on  the  inverse  of  a  given  matrix  is  in  gener¬ 
al  a  difficult  task*  It  is  hoped  that  the  present  investigation  into  the  consequences  of 
bandedness  and  total  positivity  for  the  structure  of  the  inverse  may  ultimately  prove  help¬ 
ful  in  obtaining  such  bounds* 

The  results  in  this  paper  were  obtained  in  the  study  of  a  conjecture  due  to  C.  A. 
Micchelli  [7]  •  In  connection  with  his  work  on  the  specific  approximation  scheme  of  inter¬ 
polation  at  a  (strictly  increasing)  point  sequence  X  by  elements  of  #  i*e*,  by 

splines  of  some  order  m  with  some  knot  sequence  t  -  (t^)  ,  Micchelli  became  convinced 
that  every  bounded  function  has  one  and  only  one  bounded  spline  interpolant  iff  the  partic¬ 
ular  function  which  satisfies  ffx^)  ■  (•)*  »  all  has  a  bounded  spline  interpolant  in 
$m  t  •  If  (N^)  »  denotes  the  corresponding  B-spline  basis  for  $mft  *  then 

Micchelli's  conjecture  can  be  phrased  thus:  The  matrix  A  :■  is  b°undedly  invert¬ 

ible  iff  the  linear  system 

(1)  AX  -  ((-)1) 

has  a  bounded  solution*  Micchelli  points  out  that,  for  a  finite  A  ,  this  conjecture  is 
indeed  true  and  can  be  established  using  the  total  positivity  of  A  . 

Whether  A  -  (N^(ti))  is  finite  or  not,  it  is  not  difficult  to  see  that  A  fails  to 
be  invertible  unless  A  is  m-banded,  l*e*,  unless  t  and  t  so  harmonize  that  at  roost 
m+1  consecutive  bands  of  A  are  not  Identically  zero*  It  is  then  a  small  step  to  the 
conjecture 
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A  totally  positive  m-banded  matrix  A  is  boundedly  invertible  iff  ( 1 )  has  a  bound¬ 


ed  solution" 

whose  truth  Micchelli  demonstrated  to  me  for  the  case  that  A  is  a  (biinfinite)  Toeplitz 
matrix. 

As  it  turns  out,  this  conjecture  is  incorrect;  it  fails  unless  one  assumes,  more 
strongly,  that  (1)  has  a  unique  bounded  solution.  This  uniqueness  plays  a  crucial  role  in 
the  proof  of  the  (corrected)  conjecture  given  below.  The  proof  is  first  given  for  strictly 
m-banded  matrices  and  is  then  extended  to  a  general  m-banded  matrix  by  a  limit  argument, 
using  a  'smoothing'  result  from  [4] • 

In  outline,  our  argument  is  as  follows:  We  show  that  the  nullspace 

Kl  s-  UA  s-  {f  e  R*s  Af  -  0  > 

of  a  strictly  m-banded  matrix  A  is  Haar'  in  the  sense  that,  for  every  f  eft  NO  ,  f' 
has  less  than  m  *  dlm$l  weak  sign  changes.  Here,  f*  is  obtained  from  f  by  changing 
the  sign  of  every  other  entry, 

f'U)  <-)lf<l)  ,  all  i  . 


This  makes  it  possible  to  interpolate  the  bounded  solution  x  to  ( 1 )  at  any  m-set 

I  :«  {i.,«..,i  }  by  some  yT  e  It  •  The  next  (and  hardest)  step  consists  in  showing  that, 
i  m 

for  k  -  m+  dim  with 

'il+  i“  If  et  :  lim  |f(i)|  <  »> 
i-n» 


( k ) 

yj  lies  between  x  and  0  on  the  interval  I'  9  ,  if  we  assume  that  i^  < 

...  <  .  This  implies  that  Xj.  :»  x  -  satisfies 

,xHI<k>'  «  ,x'  • 

We  use  this  fact  as  follows.  If  J  is  any  integer  interval,  and  is  the 

corresponding  section  of  A  having  the  k-th  band  as  its  main  diagonal,  then 

provided  we  choose  I  ■  [Jy(J+m)]  \(J+k)  •  But  then,  because  of  the  total  positivity  of 


A  , 


,xiij+k' 


2 


J+k  in  this  csss  .  This  uniform  boundedness  of 


and  this  is  boundsd  by  Ixl  sines  I***  ■ 
ths  invsrsss  of  all  ssetions  which  ars  principal  with  rsspset  to  the  k-th  band  is  suffic- 
isnt  for  ths  boundsd  invsrtibility  of  A  itself. 

In  this  way,  ws  show  not  only  that  sxistsncs  and  uniqusnsss  of  a  boundsd  solution  for 
(1)  implies  boundsd  invsrtibility  of  A  ,  but  gain  struct urs  information  about  ths  invsrssi 
Ths  invsrss  is  ths  pointwlss  boundsd  limit  of  ths  invsrsss  of  flnits  ssetions  principal 
with  rsspset  to  one  particular  band.  In  ths  terms  of  [1],  [2],  A  has  a  main  diagonal. 
This,  in  turn,  permits  ths  conclusior  Oat  ths  invsrss  of  a  totally  positive  band  matrix  is 
checkerboard,  a  statement  conjectured  in  [5ip.319]. 


I 


i 


3 


fV£-; 


1.  Preliminaries*  In  this  section,  we  list  certain  notational  conventions  for  easy 
reference* 

We  use  lower  case  letters  to  denote  elements  of  R1  ,  i*e*,  real  functions  on  some 
integer  set  I  ,  with  f(i)  the  value  at  i  of  the  function  (or  sequence)  f  •  If  f 
never  vanishes,  then  s(f)  denotes  the  number  of  sign  changes  in  f  ,  i.e*, 

S (  f )  |{iei:  f  (  i)  f  ( s(  i) )  <  0  }  |  , 

with  s(i)  :«  min{j  e  I  s  j  >  1}  the  successor  to  i  if  we  think  of  I  as  an  ordered 
sequence*  Here,  |J|  denotes  the  cardinality  of  the  set  J  •  If  f  vanishes  somewhere, 
then  it  is  customary  to  distinguish  between  strong  and  weak  sign  changes*  These  are  given 

by 

S  (f)  inf{S(v)  :  v  e  sign  f}  ,  S+(f)  :*  sup{S(v)  :  v  e  sign  f} 

respectively,  with  sign  f  (v  €  C-1,1}1  :  v(i)f(i)  «  |f(i)|  ,  all  i}  •  In  the  sequel, 

an  unqualified  "sign  change"  will  always  mean  "weak  sign  change".  It  is  convenient  to 

supplement  the  definition  of  s“  and  S+  by  setting 

S~(f)  :«  -1  ,  S+(f)  :«  | X |  if  f  -  0  . 

It  is  then  easy  to  check  that 

S‘(f)  +  S*(f')  -  |I|  -  1  ,  for  all  f  e  R1  , 

with  I  -  {i,,  ••*,  i}  ,  i*  <  •••  <  i_.  ,  and 
in1  n 

f’(ia)  !-  (-)8(is)  ,  . 

We  also  employ  the  prime  to  indicate  a  signature  change  in  every  other  entry  in  case 
I  *  Z  •  To  be  definite,  we  set 

f'U)  S-  (-)*£(  i)  ,  all  i  ,  all  f  e  R2  . 

In  particular,  1'  denotes  the  biinfinite  sequence  given  by 

I'd)  s-  <-)*  ,  all  i  . 

If  J  Is  a  subset  of  I  ,  then  f|j  denotes  the  restriction  of  f  e  R1  to  J  .  In 
this  connection,  \J  :«  I\J  denotes  the  complement  of  J  in  I  .  We  write  V)  instead 
of  \{j}  •  Further, 


J+k  :«  {j+k  t  j  e  J} 


0  - 


and 

[i,  j]  :=*  {k  e  Z  :  i  <  k  <  j}  . 

More  generally,  [J]  denotes  the  smallest  integer  interval  containing  J  . 

Correspondingly  (though  not  completely  in  the  same  way),  we  denote  by  AR^L  the 

restriction  of  the  "matrix"  A  6  R1  ^  to  the  subset  Kxl  of  ixj  .If  K  “  (k  #...,k  } 

1  P 

and  L  =  }  ,  with  k-  <  •••  <  k_  ,  £<...<  J l  ,  then 

Ip  1  PI  p 

det  A„  T  :=  det  ( A( k . ,£ . ) )  ^ . 

*-/*-»  13  'i,3=l 

All  norms  are  sup-norms.  Explicitly, 


Iff 


supiei | f ( i) |  ,  all  f  e  R 


and 


lAl  :*  sup{  I  Af  I  /  I  f  I  :  f  e  RJ,  lfB<«}  =  sup  E  |A(i,j)|  ,  for  A  e  R1**7 


iei  jej 


Further,  £  £  ( Z)  :=  { f  €  RZ  :  I  f  I  <  00  }  .  We  call  A  €  RZ  bounded  if  A  maps  £ 

to  itself,  or,  equivalently,  if  lAl  <  00  .If  need  be,  we  distinguish  between  the  bounded 
matrix  A  and  the  linear  map  induced  by  it  on  £w  by  calling  the  latter 

A, 


I* 


ixj 


The  matrix  A  e  R~  is  totally  positive  (  »:  tp  )  if 

det  AJC  L  >  0  for  all  KxiC  Ixj  with  |K|  -  |L|  . 

See  Karlin's  comprehensive  book  (61  for  details. 

Finally,  if  ip  and  jj,  ...,  jp  are  sequences  in  I  and  J  , 

lxJ 

respectively,  and  A  e  R  ,  then  we  use  the  customary  abbreviations 

.!«#•••# A 

jP 

and 


. £  - 


2.  The  nullspace  of  a  strictly  m-banded  bilnfinite  matrix.  The  r-th  diagonal  or  band 
of  a  matrix  A  is,  by  definition,  the  sequence  (A(i,i+r))  •  As  in  [1],  we  call  a  matrix 
A  m-banded  if  all  nonzero  entries  of  A  can  be  found  in  at  most  m+1  consecutive  bands* 
Explicitly,  the  matrix  A  is  m- banded  if 

for  some  i  ,  A(i+£,j)  ¥  0  implies  i  <  j  <  i+m  • 

Unless  otherwise  indicated  (e*g.,  by  context),  we  will  always  assume  that  £  “  0  •  For  a 
biinfinite  matrix  A  ,  this  is  merely  a  normalization  achieved  by  considering  E  A  in¬ 
stead  of  A  ,  with  E  the  shift, 

( Ef ) ( i)  :*  f ( i+1 )  ,  all  i  ,  all  f  e  RZ  , 

Z 

ar  invertible  operator  which  preserves  more  or  less  all  interesting  structures  in  R  • 

The  m-banded  matrix  A  is  called  strictly  m-banded  if 

A( i , i ) A( i , i+m)  ¥  0  ,  all  i  , 

i.e.,  the  first  and  last  nontrivial  band  is  never  zero*  In  case  of  a  biinfinite  matrix  A  , 
this  nontrivial  assumption  insures  that,  for  every  m-interval  I  «  {i+1 , • • • , i+m}  ,  every 

a  e  R1  gives  rise  to  one  and  only  one  sequence  f  with  Af  *  0  and  f  j  j  *  a  *  To  put  it 

differently,  with 

'£1:*  U  :  =  {f  e  RZ  :  Af  *  0} 

A 

denoting  the  kernel  or  nullspace  of  A  ,  strict  m-bandedness  insures  that 

for  every  I  :■  [i+1, i+m]  ,  the  map  tl  — +  Rm  :  f  I— +  is  1-1  and  onto* 

We  now  prove  this  statement  to  be  true  for  every  m-set  I  in  case  A  is  also  tp*  We 

begin  with  the  following 


then 


2.1  .  A  is  strictly  m-banded  and  tp,  and  i1  <  ...  <  ip,  j1  <  .*•  <  jp  , 


>  0 


iff  i  <  j  <  i  +m  ,  all  r  * 
-  r  r  r  - 


Proof.  Proof  of  the  *  if' 
the  customary  convention  that 


part  is  by  induction  on  p  ,  it  being  true  for  p  -  0  with 
a(J)  ■  1  *  If  i1  »  ,  then,  A  being  m-banded,  we  have 


6 


Proof*  (An  adaptation  of  the  argument  for  Theorem  5.1,2  in  [6;p.219f]).  The  assump¬ 
tion  that  p  tm  s”(xjK)  <  |J|  leads  to  a  contradiction  as  follows. 

Let  K0,  ...,  Kp  be  a  corresponding  partition  of  K  ,  i.e*,  (without  loss  of  general¬ 
ity) 

(2)  0  ?  (-)AX|K^  >  0  ,  i*0 ,  •  •  •  ,p  . 

Then 

P  i 

0  -  Z  (-)  v 

i-0 

with 

(3)  v.  Z  |x(k)|  A( • , k) ,  ,  i-0,...,p  , 

1  k«i  IJ 


7 


showing  that 
for  any  iQ  < 


is  linearly  dependent.  This  implies,  with  V  i-  [vQ 


<  i„  in  J 


V( 


Lo' 


•  •  •  #  p 


kAfflCn 
0  0 


1  |x(k  )  | ...  )x(k  )  I  A( 


V 


k  €K 
P  P 


.  |v^]  ,  that. 


•S 

,k 


) 


and  all  summands  on  the  right  are  nonnegative  by  the  tp  of  A  .  On  the  other  hand,  we  can, 
by  (2),  choose  e  ,  all  i,  so  that  |x(kQ> | • . • |x( k^) |  j*  0  and  will  therefore  have 
reached  a  contradiction  as  soon  as  we  exhibit  a  corresponding  choice  for  <  •••  <  i 
in  J  for  which 

<«»  >  0  • 

o  p 

This  we  can  do  as  follows.  Define  (i^)^1  bY 

ti-l'ip+ll  J  •  Ar  **  max^  i-i+r  '  kr“m  ^  • 

Then  iQ  <  ...  <  ip  since  both  sequences  (i^1+r)r  and  (kr-m)r  are  strictly  increasing. 

Also,  ip,  ...,  1^  €  J  since  trivially  i_^  <  ir  ,  while 

i  *  max{  i  +r  ,  k  -m  }  <  i  . 
r  -1  r  pH'1 

since  r  <  p  <  l J |  and  kr  e  K  *  [  i_«j ,  ip^+m]  •  Finally,  k^-m  <  i^  <  k^  ,  since  i_1  <  kQ 
<  ...  <  kf  implies  that  i_.j+r  <  kr  ,  hence 

kr-m  <  max{i_^+r,kr-m}  <  kf  • 

The  lemma  now  gives  (4)  and  thereby  the  desired  contradiction.  J)J 

Renark.  We  have  proved  here  a  particular  instance  of  the  statement:  "If  B  is  tp 
and  of  full  rank,  then  S  (x)  >  S+(Bx)  "  ,  provided  we  define  S*(0)  :*  length  of  0  ,  as 
we  did  earlier. 


We  conclude  that  if  y  eH  and  yjK  ^  0  ,  then 

s+(y’ |K)  -  |K|-1  -  s“(y]K)  <  |K|-1  -  |J|  -  nr-1 

hence 


Proposition  2.5 


Since  S+(f)  >  |{  i  :  f(i)  ■  0  }|  ,  this  shows  that  VI  is  then  a  Haar  space,  i.e.. 


VI  —  k"  : 

nontrivial 
else,  i.e*. 


y  I— ►  Y| ^  is  1-1  and  onto  whenever  III  =  m  •  Further,  this  shows  that  any 
y*  e  'Jt'  with  m-1  zeros  changes  sign  across  each  of  these  zeros  and  nowhere 


Corollary.  _I£  A  is  strictly  m- banded  and  tp,  and  y  e  Vl\0  vanishes  at  i,j  <  •••  < 
Vi  '  ^en  (with  iQ  :«  -»  ,  :«  ®J_<_ 

(6)  y*  (i^l  )(-)ry*  (i)  >  0  for  ir  <  i  <  ir+1  ,  r-0,...,m-1  . 

We  make  use  below  of  the  two  subspaces  xr  and  of  *  i  These  are  defined  by 

tl  {f  efl  :  lim  |f(i)|  <  •  }  . 

i+*«o 

Their  intersection  consists  of  all  bounded  solutions  to  the  homogeneous  problem  Ay  =  0  • 
The  intersection  is  therefore  trivial  iff  (0*1)  has  at  most  one  bounded  solution*  See  [2] 
for  conditions  on  A  equivalent  to  having  ft=Klf  eft'  . 
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3*  The  algorithms  L  and  R  •  Let  x  be  a  bounded  sequence  satisfying  Ax  *  1 '  •  In 


this  section,  we  investigate  y j  ,  the  sequence  in  1(1  w*1*c*1  matches  x  at  the 

m-set  X  -  {i«,..*,i  )  with  i-  <  ...  <  i_  .  For  k  -  0,**.,m,  the  interval  k  of  such  an 
m-set  is,  by  definition,  the  integer  interval 

jtk)  ^k*1  k+1*  ' 

with  i0  *-  --  ,  i^!  -  • 

Our  ultimate  goal  is  to  show  that,  for  k  ■  m+  s«  dim ,  y^*  on  1^  lies  be¬ 
tween  x'  and  0  •  The  proof  of  this  fact  involves  certain  manipulations  which  are  conven¬ 
iently  described  in  terms  of  two  algorithms,  given  below*  For  their  analysis,  the  following 
fact  is  useful* 


Proposition  3*1  •  Let  A  be  strictly  m- banded  and  tp,  and  let  s  ■  x  -  yz  for  some 
m-set  I  •  Then 

(2)  z( i)  *  0  for  all  i  0  1  . 

More  precisely, 

(3)  (-)kz'(i)  >  0  for  i  e  I(k)  ,  k»0,*..,m  . 


Proof.  Let  J  be  an  integer  interval  containing  I  in  its  ’interior'.  Since  A  is 
m- banded,  we  have 


1  •  |  J  *  ( Az )  |  j  -  Aj^L(ZjL)  ,  with  L  [J  U ( J+m) )  \  I  , 

and,  since  A  is  strictly  banded,  Aj,l  is  invertible,  by  Lemma  2*1.  Thus  (Aj^)"1  is 
checkerboard,  and 


5|l 


It  follows  that  z | L  changes  sign  strongly  at  every  slot,  with  z(j)(-)^  >  0  for  all 
j  6  J  and  to  the  left  of  I  *  Since  J  is  essentially  arbitrary,  we  conclude  that 
z'  >  0  near  -«•  and  that  Z|NI  strongly  changes  sign  at  every  slot*  ||| 


to 


Corollary*  If  the  m-set  J-{j.,.**,j)  is  obtained  from  the  m-set 

- i  m  - ' "  ■ — 

I  »  {i^,...,^}  by  moving  the  leftmost  k  &  »nts  to  the  left#  i.e.,  jr  <  ir  ' 

while  jr  -  ir  ,  r-k+1,**.,m  ,  then,  on  [i^^f  *  yj*  lies  between  x*  and  yj'  ,  with 

equality  only  at  the  points  i^+1/  •  •*,  (unless  J  -  I  ) * 

Proof*  The  assertion  follows  by  repeated  application  of  the  special  case 
jr  ■  Ar  for  r  ¥  k.  ,  ik-1  <  jk  <  ik  . 


The  proof  for  this  special  case  goes  as  follows.  By  (3), 

l-)k(yj'  -  yI')(3k)  -  (-)lc(x'  -  yI,)(jk)  <  o  . 

On  the  other  hand,  yj  and  yj  agree  at  the  m-1  points  of  ISik  ,  hence  their  differ¬ 
ence  changes  sign  strongly  across  each  point  of  ISi^  and  nowhere  else,  by  the  corollary 
to  Proposition  2*1.  Consequently, 

(4.)  <->r(yI'  -  yj’Mi)  >0  for  ir  <  i  < 

while,  again  by  (3), 

(4b)  (-)r(x*  -  yJ,)(i)  >  0  for  if  <  i  <  if+1  ,  r-k,...,m  . 

Ill 

What  is  to  follow  is  based  on  the  proposition  and  its  corollary  and  the  following 
attempt  at  constructing  a  nontrivial  element  of  * 

By  the  proposition,  S+(x* )  <  m  ,  hence  x'  has  constant  sign  near  -•  •  Let 
(5)  c  :*  sign  x*  (i)  for  i  near  -•* 

and  let  k  6  [1,m]  be  even  or  odd  depending  on  whether  e  is  1  or  -1  •  Choose  an 
n^set  X  so  that  lies  to  the  left  of  all  sign  changes  of  x*  (if  any)  . 


1 1 


Then ,  on 


1^  t  eyj*  lies  below  ex'  ,  by  (3)  •  We  then  consider  what  happens  to  y^  on  1^ 
as  we  move  . ..,  i^  to  the  left*  By  the  corollary/  eyj*  must  then  decrease.  There 

are  two  possibilities: 

(i)  No  matter  how  far  we  move  i^,  • • • ,  i^  to  the  left,  ty^*  remains  positive  on 

•  Then  we  obtain  as  a  limit  point  some  y  e  )fl  which  agrees  with  x  at  ik+1,  •••# 

,  hence  is  not  just  0  ,  and  for  which  ey'  lies  between  ex*  and  0  on  ]-»/  i^+1  [ 
But  this  implies  that  y  eH"\0  since  x  is  bounded  by  assumption.  A  refinement  of  the 
limit  process  actually  gives 

dim^T  >  m-k  . 

(ii)  Eventually,  ey^*  becomes  nonpositive  somewhere  in  1^)  f  hence  has  (at  least) 
two  sign  changes  there.  We  would  then  decrease  k  by  2  and  try  again.  By  the  corollary, 
the  two  sign  changes  of  y^ 1  just  acquired  would  not  be  affected  by  subsequent  moves. 

In  this  way,  we  either  obtain  some  nontrivial  element  of  or  else  find  ourselves 

once  again  at  (ii)  but  with  k  *  1  or  2  ,  making  further  decreases  in  k  impossible.  The 
current  y^*  must  then  have  two  sign  changes  for  every  time  we  passed  through  (ii). 

Since  S+tyj')  <  m  ,  this  limits  the  number  of  times  we  can  pass  through  (ii).  In  particul¬ 
ar,  if  we  start  with  k  *  m  or  m-1  ,  we  must  eventually  reach  (i). 

This  allows  us  to  talk  about  the  smallest  k  we  manage  to  arrive  with  at  (i)  as  we 
vary  the  initial  k  and  I  in  the  above  procedure  ?  call  it  k^  •  Analogously,  we 
define  kR  as  the  largest  k  we  manage  to  arrive  with  at  (i)  as  we  play  the  game  to  the 
right  rather  than  the  left.  The  extremality  of  kL  and  kR  and  the  fact  that  then 

dim  VC  >  m  -  kL  ,  dim  n+  >  kR 
lead  to  the  desired  conclusion  in  ways  to  be  made  precise  below. 

We  now  give  a  formal  description  of  the  game  just  played. 
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Then,  by  the  corollary  to  the  proposition, 

0  <  ey  *(i)  <  €y  .  «<i)  <  ...  <  ey  '(i)  <  ex'(i)  for  i  -r  <  i  <  i 

r  ,k  r,K+i  r  ,m  k  k 

while  yr^j(i)  *  xd)  for  1  “  ij+i  #  •••#  im  ,  and  yr^(±j)  moves  away  from  x(i^)  as 
r  increases.  Since  H.  is  finite  dimensional  and  independent  over  any  m-set  (by  Proposit¬ 
ion  2.1),  it  now  follows  that  ( yr # j ) r  has  limit  points  in  1ft.  and  any  such  limit  point 
y^  satisfies 

0  <  eyV  <  ex'  on  #  and  yj(i)  ■  X(D  *  i*ij^1  # .  •  • ,  iffl  . 

Since  x  is  bounded,  this  implies  that  zQ  yg  -  ys_f  is  in  VT  t  vanishes  at  i^, 
•«.,  i  ,  but  does  not  vanish  at  i  ,  i.e.,  the  matrix  (z  ( i  ) )  m  is  triangular 

***  o  8  t  S ,  t*  R+ 1 

with  nonzero  diagonal,  hence  invertible.  This  shows  to  be  independent  ,  hence 

dim  VC  >  m-k  •  | | | 

Finally,  if  I  is  obtained  via  EXIT3  ,  then  we  are  now  certain  that  yx  has  a  weak 
sign  change  in  each  of  the  intervals  0  ,  2,  kin  of  I  in  case  e  -  1  ,  or  in  each 

of  the  intervals  3 . kin  in  case  e  *  -1  .  This  is  so  because  once  a  sign  change  is 

obtained,  in  Step  4  ,  in  the  current  interval  k  ,  this  sign  change  persists,  by  the  corol¬ 
lary  to  the  proposition.  Further,  since  ex*  >0  on  that  interval  (we  would  have  exited 
via  EXIT1  otherwise),  it  follows  that  y^1  has  two  sign  changes  in  each  of  the  intervals 
k^n  -  2j  with  k^n  -  2j  >  0  ,  and  one  in  interval  0  if  kin  is  even,  for  a  total  of 
kin+1  sign  changes.  Since  a  nontrivial  y  e  tl  can  have  at  most  m-1  sign  changes,  an 
input  of  kin  ®  m  or  m-1  (depending  on  the  sign  of  e  )  together  with  an  Iin  which 
lies  to  the  left  of  all  sign  changes  of  x'  (to  avoid  EXIT1 )  is  guaranteed  to  bring  us  to 
EX1T2  •  In  particular,  it  makes  sense  to  define 

k^  min  {  k  s  k  obtained  as  output  via  EXIT2  from  Algorithm  L  ) 

and  then 

(6)  dim  XT  >  m  -  kL 

follows. 


Algorithm  R  is  constructed  just  as  Algorithm  L  ,  except  that  all  moves  are  made  toward 


the  right  rather  than  the  left.  For  completeness,  we  give  the  full  description* 

Algorithm  R  •  Input:  the  integer  kin  e  [0,m]  and  the  m-set  Iin  -  {ij,...,!^} 
with  i1  <  ««•  <  i^  .  Also,  set  e  :«  sign  x' (i)  for  i  near  »  • 


Step 

0. 

k 

kin  ' 

I  :> 

■  zin  * 

Step 

1* 

If 

(-)ke  < 

0  , 

or  if 

x'  changes  sign  to  the  right  of  I^k*  ,  then 

Kxrri  * 

Step 

2. 

Let 

I  :  “ 

{i.  , 

i,  +r}  ,  r-0,1,2,*..  • 

r 

1 

k 

k+1  m  *  *  • 

Step 

3. 

If, 

for  all 

r 

#  ey2  • 

>  0  on  Ir(k)  ,  then  KXIT2  . 

Step 

4* 

Pick 

an  r 

for 

which 

eyj  •  <  0  somewhere  in  Ir*k*  and  replace  I 

by  Ir  . 

Step 

5* 

If 

k  >  m-2 

,  then  BC1T3  * 

Step 

6. 

Increase  k 

by 

2  and  return  to  Step  2  • 

A  discussion  very  close  to  that  following  Algorithm  L  would  establish  the  following 
facts. 


I  —a  R  •  It.  k  is  as  on  exit  from  Algorithm  R  via  EX1T2 ,  then  dim  ft*  >  k  *  Fur¬ 
ther,  the  number 

kR  :*  max{  k  :  k  obtained  as  output  via  EXIT2  from  Algorithm  R  } 
is  well  defined,  and 

(7)  dimW  >  kR 

follows* 


we  are  now  ready  to  prove  the  main  result  of  this  section* 

Theorem  3*8  •  _I£  x  is  the  unique  bounded  solution  to  the  linear  system  Ax  *  1 '  , 
with  A  a  strictly  m-banded  bllnfinlte  tp  matrix,  then,  for  k  «  m+  :»  dim'fl*  and  for  any 
m-set  I  *  {i, •••#!}  with  ij<*»»<lm,  y^1  lies  between  0  and  x*  on  the  in- 

terval  I(k)  -  Jik,  ik+1  (  •  In  particular,  then 
(9)  |x  -  yt|  <  | x |  on  I*k)  . 
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Proof.  Let 


j0  <  •••  <  js+(x* ) 

be  points  on  which  x*  alternates  in  sign,  with  x(jQ)  ¥  0  .  Let  I  be  any  m-set  with 
ij^+r  -  jr  ,  r  -  0,1,.,., s,  and 

s  min{s+(x*),  ro-k  }  . 

L 

If  >  1  ,  then,  because  of  the  minimality  of  kL  ,  an  application  of  Algorithm  L  to  the 
input  k^-2  ,  I  is  bound  to  end  via  EXIT3,  hence  y^*  (with  a  possibly  changed  I)  has  two 
sign  changes  in  each  of  the  intervals  k^-2,  k^-4,  •  ••#  and  one  si<?n  change  in  the  inter¬ 
val  0  in  case  kL  is  even,  for  a  total  of  kL-1  sign  changes.  In  addition,  y j '  alter¬ 
nates  in  sign  on  the  points  i^  ,  ...,  i^  +8  since  it  agrees  there  with  x'  and  x' 

L  L 

does,  giving  an  additional  s  sign  changes.  We  conclude  that 

(10)  Yj'  has  at  least  1^-1  +  s  sign  changes  to  the  left  of  i^  , 

L 


and  this  conclusion  holds  trivially  in  case  kL  <  1  • 
We  now  prove  that 


Suppose  that  s  <  S  (x* )  .  Then  s  -  m-k^  *  and  we  now  know  that  has  m-1  sign 

changes  on  )-*,  i^  +g[  ,  hence  does  not  change  sign  on  [i^  +8,  *[  ,  yet  matches  x'  at 
L  L 

the  points  ik  +g,  ...,  i^  +s+(x')  on  which  x*  alternates  in  sign,  a  contradiction. 

We  conclude  that  x'  has  no  sign  changes  to  the  right  of  i^  +g  ,  hence  an  applicat- 

L 

ion  of  Algorithm  R  to  the  input  kL+s  ,  I  is  bound  to  terminate  via  EXIT2  (because  of 
(10))  with  some  k  kj  which  is  at  least  as  big  as  kL+s  ,  yet  no  bigger  than  kR  by 


the  maximality  of  kR  •  In  symbols. 


and  therefore,  with  (6)  and  (7), 


V“  *  ki  * 


dimH~  ♦  dim  >  (nt-k^)  +  k^  >  n  +  s  m  dim  +  S+(x' )  . 


This  proves  that 


dim^htL^  >  S+( x*  )  • 
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Thus,  if  now  x  is  the  only  bounded  solution  of  Ax  -  1 1  ,  then  X  (\jL  ■  10}  ,  hence 


then  s  «  S+(x')  *  0  ,  and 


dim  Tfl  >  dim  'll*  +  dim  n+ 


This  shows  that  there  must  be  equality  throughout  (12)*  This  implies  equality  in  (7),  i.e., 


and  equality  in  (11)  (with  a  *  0),  i«e., 

*L  ‘  kI  ‘  *11  ' 

and  thus  shows  (9)  to  hold  for  the  original  I  •  But  now,  since  S+(x‘ )  “  0  ,  this  could 
have  been  any  m-set  I  .  |  (  | 


Corollary  1 •  rf  A  is  strictly  oh  banded  and  t p,  and  x  is  a  bounded  sequence  satis 


fying  Ax  *  1  *  ,  then  S+(x*)  <  dim  tl“  rtUl+  • 


Corollary  2 •  The  conclusions  of  Theorem  3.8  and  Corollary  1  remain  valid  if,  in  the 
theses,  V  is  replaced  by  any  strictly  alternating  sequence  u  • 


Here,  we  call  u  atrictly  alternating  if  u(i)u(i+1)  <  0  for  all  i 


[ww  »-■»  • 


4.  The  main  result 


Theorem  4*1  .  _If  A  is  a  bounded  strictly  m- banded  biinfinite  tp  matrix,  and  x 

la  the  unique  bounded  sequence  mapped  by  A  to  1*  ,  then  A  is  boundedly  Invertible  on 

i  and  I A~ 1 1  -  Ixl  . 

00  —  — 


Proof.  For  any  integer  interval  J  ,  let  I  :*  [J  (J+m)]  (J+k)  •  Then  I  is  an 

m-set,  I  -  with  i1  <  ...  <  i^  #  say.  Let  Zj  x  -  yx  ,  with  y1  e  XL 

»  and  x  *  yx  on  I  •  Then  we  conclude  from  Theorem  3.8  that 

(2)  |zj( j) 1  <  |x( j) |  ,  all  j  e  J+k 

while 

1#|J  *  (A(x-yi))|j  *  Aj(zj|j+k) 

with 


AJ  !“  AJ,J+k  * 

Since  Aj  is  tp#  this  implies  that  Aj  is  invertible  (as  a  general  result#  though  the 
invertibility  of  Aj  could  in  the  present  circumstance  be  derived  directly  from  Lemma 
2.1),  hence  its  inverse  is  checkerboard  and  so  takes  on  its  norm  on  the  vector  1 ' i,  , 


i.e. , 


I A  I 

J 


,A/,(,V' 


J  |  J+k 


I  . 


Combine  this  with  (2)  to  get 
(3) 


I  Aj  *1  <  lxJj+kl  <  Ixl  ,  for  all  intervals  J  . 


Since  A  is  bounded  and  banded#  it  carries  c 


0 


cn ( Z)  :*  {f  e  R  :  lim 


HI- 


|f(i)  I 


o} 


to  itself#  and  the  bounded  invertibility  of  A»_  follows  now  by  a  standard  argument: 

« c0 


Let  Pj  be  the  truncation  projector, 

(P 


(f(i)  ,  i 
f)(i)  \ 

VO  oth 


€  J 
otherwise 


Then  P„ 


1  pointwise  on  cQ  ,  therefore  P  AP 


J  J+k 


A  pointwise  on  cQ  as  J 


Now  Aj  -  Aj  J+Jt  represents  the  interesting  part  of  PjAPJ+Jc  #  i.e.,  the  map  pjA|ranp 


J+k 


Therefore#  for  u  C  Cq 


and 


Uj  i*  Aj~^PjAu  e  ran  PJ+k  #  we  have 
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J+Z 


0 


,UJ  -  ■  ‘‘/’(V"  ■  pjAPj+ku)‘ 


<  I x I  IAI  lu  -  PTjLlul 
J+k 


since  IA^  I  <  Ixl  and  IP^AI  <  lAl  •  Thus  Aj^Pj  converges  point  wise 

ran  Aj^  to  a  left  inverse  of  Ai„  .  Further, 

*  c0  I co 


“1 . 


lul  -  lim  lu  I  <  lim  IA,  HAul  4  Ixl  lAul  , 

J  J 

i*e.,  A |  is  bounded  below,  hence  ran  Ai  is  closed.  The  same  argument  shows  that 


aiB°  (a|C  ) 


A  , -  is  bounded  below,  hence  ran  Ai_  is  also  dense.  We  conclude 

|c0 


0  '  1 

that  Ajc  is  1-1  and  onto  cQ  ,  hence  boundedly  invertible.  Its  inverse  is  therefore 


again  (representable  as)  a  matrix  ,  i.e.,  (A,  ) 

,c0 


-1 


•c0 


for  some  matrix  A*1  whose 


rows  are  uniformly  in  ,  and  (A|c  j"1  »  limJ>z(pjApj+k)"1pJ  pointwise  on  cQ  ,  hence 


(4) 


-1 


lijnj+z  AJ 


entrywise. 


-1  **  -1 
But  then  A  provides  the  inverse  of  A  on  i  *  ( c_  )  ,  and  I A  I  *  Ixl  since 

«  o 

-1 


IA*1!  <  Ixl  from  (3)  and  (4),  while  a"1(D  -  x  .  III 


The  assumptions  of  Theorem  4.1  can  be  weakened  in  two  ways. 

As  already  pointed  out  earlier,  the  results  of  Section  3  do  not  depend  on  having  a 
bounded  sequence  x  satisfying  Ax  *  1 1  •  It  is  sufficient  to  consider  bounded  sequences 
for  which  u  :■  Ax  is  strictly  alternating,  i.e.,  u(i)u(i+1)  <  0  ,  all  i.  For  the  results 
of  Theorem  4.1,  we  need,  more  strongly,  that  u  is  uniformly  alternating,  i.e.,  strictly 
alternating  and  with  inf  |u(i)|  >  0  •  In  that  case,  the  diagonal  matrix 

D  :  *  |  ...  ,  ( - ) Au ( i )  ,  • .  | 

is  bounded  (since  u  *  Ax  is)  and  boundedly  Invertible,  while  D^A  is  still  strictly 
m-banded  and  tp  and  carries  the  bounded  sequence  x  to  D^u  -  1'  ,  hence  D^A  is  in¬ 
vertible  on  and  |(d  *a)  *1  -  Ixl  •  Therefore  A  -  D(d”*A)  is  invertible  and 

IA*1!  <  I(D_1A)_1|Id"1I  -  «up.  . |x(l)/u(j)|  . 

I  *  J 

Secondly,  the  assumption  of  strict  m-bandedness,  though  essential  for  part  of  the  ar¬ 
gument,  is  not  essential  for  the  conclusion.  For,  according  to  (4J ,  a  bounded  m- banded  tp 
matrix  A  ,  whose  rows  and  columns  are  linearly  independent,  is  the  uniform  limit  of 
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► 


strictly  m-banded  tp  matrices  A^  (as  e  — ►  0  ,  say).  In  our  case,  the  linear 
independence  of  the  columns  follows  from  the  assumed  uniqueness  of  the  bounded  solution 
to  Ax  -  u  ,  while  the  linear  independence  of  the  rows  follows  from  the  assumed  total 
positivity  of  A  and  the  assumed  existence  of  x  with  Ax  strictly  alternating. 

Existence  and  uniqueness  of  a  bounded  solution  to  the  equation  Ax  *  u  (with  a  uniformly 
alternating  u  )  therefore  implies  existence  and  uniqueness  of  a  bounded  solution  x  to 
the  equation  A£x  *  u£  ,  with  u^  :»  u  -  (a  -  A£)x  again  uniformly  alternating  for  all 
sufficiently  small  £  •  Consequently,  A^  is  then  boundedly  invertible  on  and 

«A  -1I  <  sup  |x(i)/u  (j)|  — *  sup  I x ( i)/u( j) |  . 

£  X  ,  j  £  £^v  1  ,  J 

Thus  A  must  be  boundedly  invertible  on  £ w  and  lA  ^ I  <  sup^  |x( i)/u( j) | 

Corollary  •  The  conclusions  of  Theorem  4.1  remain  true  if  A  is  only  m-banded  and 
1 '  is  replaced  by  a  uniformly  alternating  sequence  u  • 

The  proof  of  Theorem  4.1  shows  more  than  just  the  invertibility  of  A  on  •  It 

shows  that  A  has  its  k-th  diagonal  as  main  diagonal  in  the  sense  introduced  in  [ 1 ] :  The 
sections  Aj  *  Aj  J+k  are  invertible  as  J  — ►  Z  and  the  corresponding  set  (Aj-1)  is 
bounded.  Hence  A**1  is  the  bounded  entrywise  limit  of  these  finite  matrices  Aj*1  •  Again, 
this  conclusion  persists  if  A  is  only  m-banded  since  it  is  then  the  uniform  limit  of 
strictly  m-banded  tp  matrices. 

Tlieorem  4.5  .  Let  A  be  an  m-banded  biinfinite  tp  matrix  which  is  bounded  and  bound¬ 
edly  invertible.  Then  A  has  a  main  diagonal,  i.e.,  for  some  k  and  all  intervals  J  , 

Aj  J+k  Is  invertible  and  A-1  is  the  bounded  entrywise  limit  of  (Aj  j^)*1  • 

Consequently,  with  D  the  diagonal  matrix 


<-)W1D  Is  again  tp.  In  particular,  A~*  is  checkerboard. 


(_)i+j+kA-1(i>j)  >  Q  ,  aU  j, j 


5.  Concluding  Remarks*  S.  Friedland,  in  reaction  to  a  presentation  of  these  results 


suggested  that  a  tp  matrix,  whether  banded  or  not,  must  map  onto  itself  if  its  range 

on  1 contains  1*  ,  since  it  is  then  possible  to  generate  a  pre- image  for  every  u  e 
as  a  limit  point  of  minimal  solutions  of  (PjA)y  *  PjU  ,  using  the  checkerboard  nature  of 
the  inverses  of  finite  sections  of  A  .  Further,  A.  Pinkus  showed  how  to  establish  the  sign 
regularity  of  DA-1D  ,  with  A-1  the  bounded  inverse  of  a  tp  matrix  A  ,  without  assuming 
that  the  inverse  is  the  limit  of  inverses  of  finite  sections*  These  matters  are  made 
precise  in  [3]* 
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